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Abstract 

We perform molecular dynamics simulations to investigate the diffusive motion of TIP5P water 
in the direction perpendicular to the two hydrophobic confining walls. To calculate the diffusion 
coefficient, we use the concept of the characteristic residence time which is calculated from the 
exponential decay of the residence time probability distribution function. We find that a diffusion 
anomaly of water, increase of diffusion upon compression, is absent in the direction perpendicular to 
the confining walls down to the lowest temperature we simulate, 220K, whereas there is a diffusion 
anomaly, similar to that in bulk water, in the direction parallel to the walls. The absence of a 
diffusion anomaly may arise mainly due to nanoconfinement, rather than due to the hydrophobic 
property of the confining walls. 

PACS numbers: 66.10.C-, 61.20.Ja, 66.10.-x 
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In addition to the thermodynamic anomahes of water 
dynamic anomahes. For example, in contrast to simple liquids, where the diffusion coeffi- 
cient D decreases upon compression, for water at low temperature D first increases upon 
coinpression until reaching a maximum and then decreases upon further increase of pressure 

Confined water shows many structural and dynamic properties different from bulk water 
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confinement and t he interaction with the confining walls 

261. Studies of water confined 
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between two hydrophobic plates suggest that the diffusion coefficient D\\ along the direction 
parallel to the plates displays a diffusion anomaly, which is the existence of a maximum of 
D\\ as a function of a density at constant temperature, similar to bulk water i23|. Moreover, 
the temperature below which the anomalous diffusive behavior occurs is shifted to lower 
temperature by about 40 K compared to bulk water 2^. 

Dynamics in the direction perpendicular to the confining plates is modified by the small 
distance between plates (~ 1 nm). Studies have been carried out to investigate the diffusive 



motion in the confining direction for many confined systems [27|, l28|, l29|, l30|] . Liu et al. showed 
that a different treatment for the diffusion coefficient is needed in the confining direction 
and used a "dual simulation method" to calculate the more precise diffusion coefficient at 
a liquid- vapor interface {2^. Sega et al. pointed out two different time regimes (100 ps 
and 1 ns) of the mean square displacements (MSD) for the complete description of diffusion 
in confined water and showed a nonlinear MSD due to both spatial inhomogeneity and 
confinement The diffusion coefficient D_l in the direction perpendicular to the walls 

is difficult to find precisely in nanoconfined water due to the finite and very small space 
available for water molecules. Before reaching the diffusive regime at which MSD is linearly 
proportional to time, MSD already enters the plateau regime due to displacement bounded 
by a finite nanosize space. Hence, D_l cannot be extracted from the MSD. The Green-Kubo 
relation for D± is also not valid for the confined system, as mentioned in Ref. 271]. 

Here we propose an alternative approach to calculate and ask whether the diffusion 
anomaly of water, which has been shown to exist in the parallel direction and in bulk 



dynamics, also exists in the direction perpendicular to the hydrophobic confining 



We perform molecular dynamics (MD) simulations of = 512 TIP5P 3l|] water 



ates. 



molecules confined between two infinite parallel smooth plates, separated by 1.1 nm, which 
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are able to contain 2~3 layers of water molecules. The plates are located at the positions z= 
±0.55 nm. We model the water-wall interaction by a 9-3 Lennard- Jones (LJ) potential which 



IS common 
plates 
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y used to represent the effective interaction of water molecules with the confining 



33] . We choose the parameters for the confining potential to have the hydropho 



bic property, as in Ref. [23i]. Details of the simulations are also described in Ref. [23|]. We 
perform MD simulations for seven temperatures T = 220, 230, 240, 250, 260, 280, and 300 K 
and for eight densities p = 0.80, 0.88, 0.95, 1.02, 1.10, 1.17, 1.25, and 1.32 g/cm^. These 
densities are calculated by considering the accessible space between the walls, as explained 



in Ref. [23|]. Periodic boundary conditions are used in the x and y directions. 

To calculate D_i_, we first divide the system into three residence regions along z direction 
such that there are two symmetric adjacent regions to the surface and one middle region 
(see Fig. [1]). The width of each region is 0.14 nm, and the separation between two adjacent 
regions is Rz = 0.28 nm which is the same as the linear size of a water molecule. Then, we 
calculate the residence time distributions P{tr) of water molecules in the given residence 
region. The residence time is defined as the time over which water molecules stay in one 
region before leaving it. 

In Fig. [21 we show that P{tr) of the hydrophobic confined water decays exponentially for 
all temperature investigated. It was shown that the anomalous dynamics of water becomes 



subdiffusive when it is confined in Vycor pores with the condition of low hydration [3J]. This 
subdiffusive motion is related to a power law behavior of P(rij). However, our results show 
that in hydrophobic confinement, P{tr) has an exponential behavior for all temperatures 
and densities investigated. We calculate the characteristic residence time t'^ by finding the 
inverse slope of a straight line fit to P{tr) on the semi-log plot 

PirR) - exp . (1) 

On average, during t'j^ water molecules diffuse the same distance as the separation of two 
defined regions, as shown in Fig. [H 

Since water molecules, on average, diffuse the distance equal to the separation between 
two regions in the time interval t'^, we can write down the diffusion coefficient D± 

where D± is the diffusion coefficient in the direction perpendicular to the walls, (Rz) denotes 
the separation between two residence regions, and {t^) is the characteristic residence time 



= (2) 



averaged over an ensemble. To get the average value of Dj,, we use the characteristic 
residence times averaged over three different residence regions. As shown in Fig. [2](b), 
P(r/j) has the same exponentially decaying behavior in different residence regions for all 
temperatures and densities investigated. 

Since the calculation of D\\ is possible from the Einstein relation |3a], we investigate the 
validity of Eq. (2) by calculating D\\ using both the Einstein relation and the characteristic 
residence time. In Fig. |3], we show D\\ as a function of density for T=240 K calculated from 
both methods for a comparison. We find that the method using the characteristic residence 
time gives slightly larger value of D\\ than using MSD, as shown in Fig. [3l Both methods 
for D\\ in our simulations exhibit a diffusion anomaly with a maximum at p = 1.02 gjcw? 
as a function of density, same as in Ref. 23|]. Therefore, we can say that Eq. (2) gives the 



correct value of D and can be used for the calculation of D_l to investigate the existence of 
a diffusion anomaly in the perpendicular direction. 

In Fig. m we show D ^ as a function of density for all temperatures studied. D ^ decreases 
as the density increases. Contrary to the diffusion anomaly found in the parallel direction 
(similar to bulk water), as shown in Fig. [3], our results show a diffusion anomaly in the 
perpendicular direction does not exist down to the lowest temperature we simulated. As a 
result, we conclude that a diffusion anomaly of water is absent in the confining direction down 
to very low temperatures in nano confinement. From the fact that a diffusion anomaly exists 
in the parallel direction but not in the perpendicular direction, the main contribution to the 
absence of a diffusion anomaly in the perpendicular direction might be the nano confinement 
rather than the hydrophobic property of the confining walls. 

Next we study D ^ as a function of temperature along a constant density path. In Fig. [5l 
we find the temperature dependence of can be fit with a Vogel-Fulcher-Tammann (VFT) 
form for all the densities studied, 

A 



= L)° exp 



(3) 



T-To, 

Here D^, A, and Tq are fitting parameters, and we use 160 K as the value of the parameter 
Tq. It has been experimentally observed that there is a fragile-to-strong transition near T= 
220 K in both supercooled water confined in micellar templated mesoporous silica matrices 



MCM-41 



36| and DNA and protein hydration wate r |37|, w 



to a liquid-liquid phase transition scenario 
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lich is shown to be connected 
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431]. In computer simula- 



tions of TIP5P hydration water at atmospheric presseure, this non-Arrhenius to Arrhenius 
crossover occurs at T ^ 250 K, the temperature at which the isobaric specific heat has a 
maximum 
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42]. In contrast to simulation results of T1P5P bulk water, we find that 



in the confining direction the TIP5P water confined between two hydrophobic plates does 
not show a non-Arrhenius to Arrhenius dynamic crossover down to the lowest temperature 
we simulated, T = 220 K, suggesting that if there is a crossover in the dynamics it would 
occur at much lower temperature compared to bulk water. This finding is consistent with 
the temperature shift found for thermodynamic and dynamic properties of water confined 
between hydrophobic surfaces |23l. 1441]. 

In summary, we have performed MD simulations of TIP5P water to calculate and investi- 
gate the self-diffusion coefficient in the direction perpendicular to the hydrophobic confining 
walls. We calculated the diffusion coefficient using a form similar to the Einstein relation 
of self- diffusion, described with the separation between two residence regions and the char- 
acteristic residence time found from the exponentially decaying residence time distribution. 
Contrary to the diffusive dynamics of bulk water and the dynamics of confined water in the 
parallel direction, our simulation results show that a diffusion anomaly does not exist in 
the direction perpendicular to the confining walls. As the density increases, the diffusion 
coefficient keeps decreasing without reaching a maximum over the whole temperature range 
investigated. By comparing the perpendicular direction to the parallel direction, it is reason- 
able to conclude that this absence of a diffusion anomaly may arise due to nano confinement 
rather than the hydrophobic property of the confining walls. In addition, we find that the 
temperature dependence of the diffusion coefficient along the constant density path in the 
perpendicular direction shows a VFT form. 
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grant CIIE0606489 for support. 
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FIG. 1: (color online). Schematic description of the motion of water molecules in confined space 
with three defined residence regions. Two confining walls are located at z = ib0.55nm. The red 
circle represents a water molecule. We define the size of one residence region as 0.14 nm and 
the separation between two regions as 0.28 nm. t'^ denotes the characteristic residence time 
calculated from the residence time distribution. On average, water molecules diffuse the distance 
of the separation between regions, Rz = 0.28 nm, in the z direction perpendicular to the confining 
walls during a time t^. 
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FIG. 2: (color online), (a) A semi- log plot of the residence time distribution function P{tr) of 
TIP5P water molecules at several temperatures and density p = 0.80 g/cm^. For all temperatures 
and densities investigated here, P{tr) decays exponentially [~ exp(— rij/r^j*)]. (b) The residence 
time distributions P{tr) for different residence regions (see Fig. [TJ) at /) = 0.80 g/cva? and 
T = 250 K and 300 K. 
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FIG. 3: (color online). The diffusion coefficient D\\ in the direction parallel to the confining walls 
as a function of density. D\\ is calculated by using both from the Einstein relation (MSD) and the 
charateristic residence time (RT) at T=240 K. Both calculations exhibit a diffusion anomaly in D\\ 
as a function of density with a maximum at p = 1.02 g/cm? . 
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FIG. 4: (color online). The diffusion coefficient D± in the perpendicular direction as a function 
of density along constant temperature paths. Z)j_ decreases as density increases over the entire 
temperature range investigated. Our results suggest that there is no diffusion anomaly along the 
perpendicular direction even at the lowest temperature T=220 K we simulated. In the parallel 
direction, however, there is a diffusion anomaly at temperature lower than T = 250 K. (see Fig. [3] 
and Ref. I23].) 
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FIG. 5: (color online). A semi- log plot of -D_l as a function of l/(r — Tq) along the constant density 
paths. The temperature dependence of D±^ can be fit with a Vogel-Fulcher-Tammann (VFT) form 
using the parameter value of Tq = 160 K. 
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